High- Tc ferrolectricity emerging from magnetic degeneracy in cupric oxide 
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Cupric oxide is multiferroic at unusually high temperatures. From density functional calculations 
we find that the low-T magnetic phase is paraelectric and the higher-T one ferroelectric, with a size 
and direction of polarization in good agreement with experiment. By mapping the ah initio results 
onto an effective spin model we show that the system has a manifold of almost degenerate ground 
states. In the high-T magnetic state non-coUinearity and inversion symmetry breaking stabilize each 
other via the Dzyaloshinskii-Moriya interaction. This leads to a novel mechanism for multiferroicity, 
with the particular property that non-magnetic impurities enhance the effect. 

PACS numbers: 71.45.Gm, 71.10.Ca, 71.10.-w, 73.21.-b 



In multiferroics the simultaneous presence of electric 
and magnetic ordering is particularly intriguing when the 
magnetic ordering triggers the ferroelectric polarization, 
as was observed for the first time by Kimura and cowork- 
ers in TbMnOs [l[. Since then, several so-called type- 
II multiferroics |2l| have been discovered in which mag- 
netic order causes ferroelectric order. Although plenty 
of potential applications are envisioned, in random ac- 
cess memory devices for instance, the small values of the 
induced polarization as well as a low transition temper- 
ature in most type-II multiferroics hinder practical ap- 
plications. The very recent discovery that cupric oxide 
(CuO) is a type-II multiferroic with a high antiferromag- 
netic transition temperature T^ of 230 K changed this 
situation drastically and opened the perspective to room- 
temperature multiferroicity [Sl \^ . The discovery is even 
more intriguing considering that CuO is closely related 
to the family of copper-oxide based materials displaying 
High- Tc superconductivity. 

From a theoretical point of view, the microscopic mech- 
anism of multiferroicity in CuO is not clear yet, partic- 
ularly because its type-II behavior is apparently not a 
groundstate property: it is only present at finite temper- 
atures, between ~210 and 230 K, disappearing above and 
below. Here we clarify the mechanism for the observed 
finite temperature multiferroicity in CuO. 

To elucidate this point we have investigated the elec- 
tronic structure of CuO with density functional calcula- 
tions for the different magnetically ordered phases. These 
calculations confirm, as we will see, the presence of mag- 
netically induced ferroelectricity in CuO and we find 
a polarization that agrees with experiment. A subse- 
quent investigation of the stability of magnetic phases 
at finite temperatures using classical Monte-Carlo sim- 
ulations shows that the experimental ground state at 
low temperature can be well understood by mapping 
the magnetic interactions onto a Heisenberg Hamiltonian 
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FIG. 1: (Color online) Schematic view of API and AF2 mag- 
netic states. Cu ions belong to different plane along y = G 
(red/big circle, termed even plane in the text) and y = 1/2 
(green/small circle, termed odd plane in the text). Dot (cross) 
refers to spin pointing along the positive (negative) y axis. 



with ah initio derived exchange constants. The micro- 
scopic model shows that the ground state is almost de- 
generate so phase selection occurs through small terms. 
We show that multiferroicity in CuO arises from a new 
mechanism in which spin canting and polarization mutu- 
ally stabilize each other, crucially involving phase selec- 
tion through the Dzyaloshinskii-Moriya interaction. 

We have studied the electronic structure of CuO by 
performing calculations using the PAW method as imple- 
mented in VASP [5[. To take into account the Coulomb 
interactions between the Cu id electrons we employ 
SGGA-t-U W^ and hybrid functional (HSE) 1] schemes 
for Ueff = U — Jh ranging between 3.5 and 7.5 eV, 
and fraction of Hartree-Fock (HF) exchange (a) between 
0.15 and 0.25. In the experimental structure of CuO 
(C2/c space group No. 15), Cu ions are arranged as cor- 
ner and edge-sharing square-planar CUO4S, in which Cu 
ions lie in a plane formed by oxygen neighbors forming 
0-centered tetrahedra [10|. The low temperature mag- 
netic structure (AFl) with Cu magnetic moments aligned 
collinearly along the y axis and ordered antiferromagnet- 
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TABLE I: Exchange coupling parameters (meV) calculated 
within SGGA-I-U and hybrid functional calculations. The 
structure allows for Ja 7^ Jd and Jb 7^ Jc but we take them 
equal for simplicity. This is inessential for our conclusions. 
We keep the same notation of reference [la ]. 



ically along z and ferromagnetically along the x direction 
(see Fig. [T|) is found to be the ground state in agreement 
with the neutron diffraction study 111 and other ab-inito 
calculations [18[ |l9|. The high Neel temperature sug - 



gests the presence of strong exchange interactions [12 



14[ | and the quasi-one-dimensional antiferromagnetism of 



CuO is suggested by experiments [15|, [16|. The energy 
gap of 1.4 eV [l7| is recovered using Ue£f=5.5 eV or an 
exact- exchange fraction as a = 0.15 [l8|. The Cu and 
O spin moments are 0.6 /is and 0.1 fj,B respectively, in 
good agreement with the results of neutron diffraction 
experiments [lli]. 

We analyze the exchange interactions using the Heisen- 
berg Hamiltonian Hm = J2ij Jij ^i ' ^j- The parameters 
Jij , reported in Table [H are calculated from total energy 
differences of different magnetic configurations within a 
unit cell containing 32 Cu sites. The strongest interac- 
tion is Jz ~ 100 meV in good agreement with the optical 
experiments of Ref. [iGj and the neutron data of Ref. [ij] . 
The magnetic structure within a constant-y plane can be 
viewed as chains running along z with dominant inter- 
action Jz and moderate interchain interactions Jx, J2a 
and J26. Hereafter we will term the planes with 2y odd 
(even) as odd (even) planes. Jy couples planes of the 
same kind (separated by Ay = ±1). Finally J^, J;,, Jc 
and Jd couple nearest neighbor (nn) planes of different 
kind. 

Remarkably the classical coupling energy among un- 
like planes vanishes when the chains are assumed to be 
aligned AF. Indeed Jy favors a ferromagnetic staking of 
planes of the same kind. A close examination of the struc- 
ture shows that the pattern of couplings among the plane 
at y = 1/2 and the plane at y = 1 is identical to the one 
of Fig. [1] but with the exchanges J a o Jd and Jc <-> Jb- 
It is easy to check that this symmetry makes the classical 
energy of the model independent of the angle between the 
magnetization in even and odd planes, and the ground 
state is infinitely degenerate. This degeneracy will play 
an important role for the multiferroic mechanism. 

Relaxation of charge and lattice and small anisotropies 
favor a collinear structure. Indeed within DFT calcula- 
tions, including spin-orbit (SO) coupling, we find that 
AFl is the ground state with an energy gain of 2.2 nieV 
per Cu atom respect to the non-coUinear AF2 structure 
shown in Fig. [TJ The easy axes is found to be along the 



y direction. 

In CuO, an incommensurate phase with magnetic 
modulation vector Q=(0.506, 0, 0.517) has been reported 
at the temperatures where the ferroelectricity is ob- 
served [3I, [ll|, [l2|. In this phase spins order along the 
y axis and in the xz plane with spins on even and 
odd planes perpendicular to each other. This is shown 
schematically in the right panel of Fig. [T] where the small 
incommensuration of the structure has been neglected. 
Experimentally, a small electric polarization P (^0.01 
[iCjcrr?) is found along y axis[3|. 

On a first sight one would expect that the incommen- 
surate spiral is crucial to obtain a finite polarization, as 
in the standard cycloid scenario [j, [20|, [21 1. However the 
commensurate state closest to the incommensurate spi- 
ral, labeled AF2 in Fig. [1] has spin canting which can 
produce a finite polarization. Indeed, taking into account 
spin-orbit coupling, we evaluate the electronic contribu- 
tion to the polarization P using the Berry phase (BP) 
method [22| on the commensurate AF2 state and we ob- 
tain Paf2 ^ 0.02 fiC/cm^ along y axis, in overall good 
agreement with the experimental value. Thus the incom- 
mensurate state is not crucial but canting clearly is. The 
perpendicular configuration ensures that AF2 state has 
maximal spin current ji_2 = (Si x S2) among nn planes 
of different kind. We will show that this is a fingerprint 
of the proposed scenario. 

While canting and spin orbit coupling are standard 
ingredients of the cycloid scenario causing the multifer- 
roicity [J, [20|, [21| , the situation in CuO is subtly different. 
In order to explain the difference and similarities we il- 
lustrate the two mechanisms in the one-dimensional (ID) 
model [20| depicted in Fig. [5] Consider an hypothetical 
Cu-0 chain with Hamiltonian H = Hm + Hdm + He- 
For the purely magnetic part Hm we assume there is a 
nearest neighbor AF interaction Ji and a next nearest 
neighbor AF interaction J2- The Dzyaloshinskii-Moriya 
(DM) interaction [2J, [25| and elastic contributions read: 

Hdu — 2^ A(u„_|_i/2 >< e„^„+i).(S„ X S„+i), (1) 



He - 2_/ o'^^+i 



/2I 



Here u are the oxygen displacements and e„_„_|_i is a unit 
vector joining nearest neighbors atoms. Treating the spin 
classically for J2 > Ji/4 one finds that the ground state 
is a spiral with a pitch angle given by cos 61 = —J1/AJ2 
and a finite spin current j = (S„ x S„+i). SO interac- 
tion is not necessary to stabilize this state [Fig.[5Ja)]. In 
the presence of SO coupling the free energy per site due 
to uniform displacements u of oxygens, perpendicular to 
the chain and to j is given by SFdm = ^uj + fu^- Min- 
imizing, one obtains a polarization P — Squ — ~Sqj/k 
[Fig. [2Jb)] due to the difference in charge Sq of the two 
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FIG. 2: (Color online) (a) Schematic mechanism of multifer- 
roic effect in ID. a) and b) standard cycloid scenario as in 
Ref. [2y|. In c) Ji = and the system separates into two in- 
terpenetrating sublattices. The angle between the magnetic 
moments of the different sublattices is undefined, d) SO in- 
teraction stabilizes both the angle and the finite polarization. 



In the case Ji — [(c) and (d)] the system separates 
into two interpenetrating sublattices analogous to the 
odd and even planes of CuO. Without SO coupling the 
angle between the magnetic moments of the different sub- 
lattices is arbitrary and the ground state is infinitely de- 
generate (c) . This degeneracy can be broken by the DM 
interaction. The free energy per site can be expanded as: 
SF = \xjj'j'^ + -^""J + 1"^- Here Xjj is a spin current 
susceptibility defined for Hm alone. 5F has to be min- 
imized with respect to both j and u since Hm does not 
determine the spin current. Minimizing with respect to 
u one obtains 6F = Hxjl - ^^ /k)p. When 
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it is convenient to maximize the spin current j oc sin0 
where (j) is the angle between the magnetization in dif- 
ferent sublattices. Thus the energy acquires a term 
5F oc — sin which favors perpendicular magnetic mo- 
ments on nearest neighbor sites [Fig. [2ld)]. At the same 
time this gives rise to a finite polarization as for the cy- 
cloid mechanism. Notice however that the spontaneous 
breaking of symmetry and the spin canting drive each 
other unlike the cycloid scenario where the spin canting 
is driven by the magnetic Hamiltonian. From Eq. ([2]) we 
see that the multiferroic effect is favored by strong SO 
coupling, soft lattices and a large spin current suscepti- 
bility. This effect competes with thermal and quantum 
fluctuation which favor a collinear configuration, accord- 
ing to the order by disorder mechanism, and tend to sup- 
press Xiil23|. Also coupling to other lattice distortions 
will favor a collinear state trough effective biquadratic 
terms in the Hamiltonian [~ (S^.S^)^] and suppress Xjj- 

The ID mechanism can be easily generalized to CuO 
by replacing each magnetic site in Fig. [5] by a constant-y 
plane of Cu atoms shown in Fig. [TJ 

In order to estimate Xjj for the CuO structures we per- 
form classical Monte-Carlo simulations on the following 
3D spin Hamiltonian: H = Hm + Tliij) ^ ' J«j- Here, 
the first term is the magnetic Hamiltonian of CuO and 



the second term describes the linear coupling of the spin 
current among unlike planes to an auxiliary external field 
D. Summation index (ij) represents inter-plane nearest 
neighbor bonds indicated as J a, Jb, Jc and Jd in Fig. 
[H We take D directed along x-axis for all pairs in the 
unit cell. Rather than measuring the spin current we 
compute the sum of the relevant components which is 
proportional to the polarization, p — J2^ij ^ i^i ^ ^j)- 
The desired susceptibility is given by Xjj — Py/Dx in the 
limit of vanishing D^ ■ When Xjj is large and the condi- 
tion Eq. ^ is satisfied a spontaneous D and polarization 
will stabilize each other as explained above. Notice that 
D breaks inversion symmetry in the system whereas the 
high temperature structure of CuO has inversion so there 
is no "permanent" D in the high temperature phase. Of 
course there will be DM couplings at high-temperatures 
in the CuO structure but those preserve inversion sym- 
metry and are not related to the appearance of the polar- 
ization. For simplicity we neglect these latter couplings. 

To lift the degeneracy between AFl and AF2 in favor of 
the former we introduce a weak anisotropy in the Heisen- 
berg exchange term. This is done by replacing J^S^ • Sj 
by J.iSfS^ + (1 + j)SfS] + S^S^) in Hm- Other terms 
like the biquadratic contribution are expected to have a 
similar effect. We use 7 — 0.02 which translates into an 
anisotropy energy of ~ 7 Jz = 2.15 meV. 

We employ a classical Monte Carlo (MC) technique 
to explore the competition between different magnetic 
states at zero and finite temperatures. Given that the 
CuO is a system with spin 1/2, the quantum effects in 
this system are unavoidable. Nevertheless the interesting 
transitions occur at high temperatures where it is safe to 
assume a classical renormalized regime [26J. In order to 
simplify the Monte Carlo computation we consider only 
4 possible states at 90° for the spin variables. Figure 
[3l^a) shows the phased diagram in the T-D plane. In the 
absence of the external field D, the system undergoes a 
transition from a paramagnetic (PM) to an AFl state 
with T/v ^ 2bQK . Presence of a small D opens a narrow 
window near T^ where AF2 is stabilized. A large exter- 
nal field eventually drives the groundststate to be AF2 
for D = 7Jz/8 - 0.27 meV . The inset of Fig.|3i;b) shows 
the susceptibility computed as the ratio Py/D^ for small 
Dx- Remarkably a strong peak appears around the PM 
to AFl transition which will favor a spontaneous polar- 
ization of the system. 

In order to check the mechanism we again consider the 
3D CuO model with the terms Hdm and He analogous 
to the ID model. Assuming classical lattice displace- 
ments they can be integrated out of the partition func- 
tion leading to a quadratic effective interaction among 
spin currents so H^m + He is replaced by Hdme — 
-(AV2fc)E(jj)(S» X Sj)2. Fig. Eljc) shows the sponta- 
neous polarization p as a function of temperature for the 
model defined by H — Hj^j + Home- As expected one 
finds that, close to the PM to AFl transition, Hume 
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FIG. 3: (Color online) (a) Phase diagram of the magnetic 
model of CuO with an external field coupling linearly with 
the spin current, (b) Polarization of CuO model with the ad- 
dition of a biquadratic spin-current term obtained eliminating 
the lattice degrees of freedom in the DM coupling. Curves are 
labeled by the value of X^ /{2k). The inset shows the suscepti- 
bility computed as Py/D. Within linear response one should 
take the limit D — )• corresponding to the upper curves. 



stants derived by ab initio simulations we also confirm 
the high T/v of this compound. We explain the multifer- 
roic efi'ect as arising from a new mechanism in which spin 
canting and polarization mutually stabilize each other 
with a crucial role of Dzyaloshinskii-Moriya interaction. 
Our results open new routes for the material design of 
multiferroics. 
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203523) and IIT-Seed project NEWDFESCM. We thank 
A. Boothroyd for a critical reading of the manuscript. 

During the completion of this manuscript we became 
aware of Ref. j28|| , presenting results from ab initio calcu- 
lations similar to ours. The mechanism for multiferroicity 
that we present here, however, is very different from the 
findings in Ref. 28]. 



induces a phase with broken inversion symmetry and a 
spontaneous polarization as seen in the experiment. The 
peak in the polarization is very similar to the experimen- 
tal observations of a finite electrical polarization between 
230K and 213K ^3]. If the parameter X"^ /{2k) is made too 
large (> 1 meV) the ferroelectric phase extends to zero 
temperature. 

According to our Monte Carlo simulations the main 
reason for the spontaneous polarization is a strongly 
enhanced spin-current susceptibility close to the AFl- 
PM transition. While there is of course a divergent 
staggered susceptibility when approaching the AFl-PM 
phase transition the enhancement of the unrelated spin- 
current susceptibility is not trivial. This is similar to 
the physics of quantum critical points relevant to heavy 
fermion compounds 27| where close to the quantum tran- 
sition between a disordered an a magnetically order state 
a different order appears (superconductivity) which can 
be attributed to an enhanced pairing susceptibility. 

As mentioned above thermal and quantum fluctua- 
tions tend to suppress Xjj ^^^d the polarization. On the 
other hand disorder on the magnitude of the magnetic 
moments will enhance the tendency to have perpendic- 
ular orientations among the sublattices and enhance the 
polarization |23| which opens a new way to engineer high- 
Tc multiferroic materials. Non-magnetic impurities, va- 
cancies or magnetic impurities with a different spin will 
lead to this effect. In addition our results suggest to 
search for other materials where two subsystems have 
negligible interactions by symmetry but strong interac- 
tions within one subsystem as a recipe to discover new 
manipulable multiferroics. 

To summarize, our density functional calculations con- 
firm the magnetically induced ferroelectricity in CuO 
with polarization in agreement with experiments. By 
combining Monte-Carlo analysis with the exchange con- 
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